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Carbon  nanotubes  have  been  identified  as  promising  reinforcing  materials  for  high- 
performance  nanocomposites  (e.g.,  Thostenson  et  al.,  2001;  Maruyama  and  Alam,  2002). 
Reliable  structural  applications  of  carbon  nanotube-reinforced  composites  depend  on 
accurate  understanding  of  their  mechanical  behavior.  Efforts  have  been  made  to 
characterize  the  mechanical  properties  of  the  said  nanocomposites.  These  studies,  being 
typically  based  on  experimental  measurements  or  molecular  dynamics  simulations,  tend 
to  be  expensive  and  configuration/material  specific.  The  use  of  continuum-based  models 
can  mitigate  these  difficulties  (e.g.,  Odegard  et  al.,  2002;  Li  and  Chou,  2003a;  Gao  and 
Li,  2003)  and  is,  therefore,  very  desirable.  The  goal  of  the  current  research  project  is  to 
explore  continuum-based  models  for  predicting  effective  elastic  properties  of  and 
describing  interfacial  load  transfer  in  carbon  nanotube-reinforced  polymer  composites. 

A  continuum-based  micromechanics  model  has  been  proposed  (see  Li  et  al.  (2004)  for 
details).  Two  steps  are  involved  in  the  modeling  procedure  for  two  different  length 
scales.  In  the  first  step,  which  is  at  the  nanometer  scale,  a  single  walled  carbon  nanotube 
(SWCNT),  the  non-bulk  polymer  adjacent  to  the  nanotube,  and  the  nanotube/polymer 
interfacial  layer  are  collectively  treated  as  an  effective  continuum  fiber  (isotropic  or 
transversely  isotropic),  whose  effective  elastic  properties  are  determined  by  using  an 
equivalent-continuum  modeling  method  (Odegard  et  al.,  2002).  In  the  second  step,  which 
is  at  the  micron  scale,  the  Mori-Tanaka  method  restructured  by  Tandon  and  Weng  (1984, 
1986)  and  Qiu  and  Weng  (1990)  is  used  to  predict  the  effective  elastic  properties  of  the 
SWCNT-reinforced  polymer  composites.  As  a  parametric  study,  four  nanocomposite 
systems,  two  transversely  isotropic  ones  containing  unidirectionally  aligned  nanotubes 
and  two  isotropic  ones  having  randomly  oriented  nanotubes,  are  considered.  Nanotube 
aspect  ratio  (AR),  nanotube  volume  fraction  (ci)  and  nanotube  orientation  are  chosen  as 
three  controlling  parameters.  Thermoplastic  polyimide  LaRc-SI  with  a  3%  molecular 
weight  offset  is  taken  as  the  isotropic  matrix  material.  Its  Young’s  modulus  and  Poisson’s 
ratio  are,  respectively,  Eo  =  3.8  GPa  and  vo  =  0.4  (Whitley  et  al.,  2000).  For  the 
transversely  isotropic  effective  fiber  the  five  independent  elastic  stiffness  constants  are, 
respectively,  taken  to  be  Cn  =  457.6  GPa,  Cn  =  8.4  GPa,  C22  =  14.3  GPa,  C23  =  5.5  GPa, 
C44  =  27.0  GPa  (Odegard  et  al.,  2003).  For  the  isotropic  effective  fiber  the  Young’s 
modulus  (E\)  and  Poisson’s  ratio  (vi)  are,  respectively,  regarded  to  be  the  same  as  the 
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longitudinal  Young’s  modulus  and  Poisson’s  ratio  of  the  transversely  isotropic  effective 
fiber,  i.e.,  E\  =  E{[  =  450.4  GPa,  V|  =  v{2  -  0.42,  where  EJU  and  v{2  can  be  obtained 
from  the  stiffness  constants  given  above.  The  effective  fiber  is  assumed  to  have  the  same 
length  as  that  of  the  nanotube,  and  the  nanotube  volume  fraction  (c i)  is  taken  to  be  34% 
of  the  effective  fiber  volume  fraction  (Odegard  et  ah,  2003). 

For  the  transversely  isotropic  composites  containing  unidirectionally  aligned  isotropic 
nanotubes,  numerical  results  show  that  when  AR  increases  the  effective  longitudinal 
Young’s  modulus  E\\  increases,  while  the  effective  Poisson’s  ratio  vj2  decreases,  both  in 
a  monotonic  manner  (see  Figs.  1  and  2).  Furthermore,  with  the  increase  of  AR,  £11 
approaches,  while  v\2  deviates  from,  their  corresponding  results  predicted  by  the  rule  of 
mixtures  (e.g.,  Chawla,  1993)  according  to: 

£11=£0+(£l-£0)c,/0.34,  (1) 

vl2  =  vo+(vl-vo)cl/034.  (2) 

In  addition,  v\2  increases  with  c\  when  AR  is  small,  but  decreases  with  c\  when  AR  is 
getting  close  to  25.  When  isotropic  nanotubes  are  randomly  oriented  in  the  polymer 
matrix,  the  effective  Young’s  modulus  E  of  the  resulting  isotropic  composites  increases 
with  both  AR  and  c\,  while  the  opposite  trend  is  observed  for  the  effective  Poisson’s  ratio 
v(see  Figs.  3  and  4).  By  comparing  Fig.  1  with  Fig.  3  and  Fig.  2  with  Fig.  4,  one  can  see 
that  the  values  of  E  and  vare  smaller  than  their  corresponding  values  of  Eu  and  vn  for 
the  same  values  of  AR  and  c\,  as  expected. 


Vo  1  urn*  fraction,  cl{%) 


Fig.  1.  Aligned  SWCNTs  (isotropic).  Fig.  2.  Aligned  SWCNTs  (isotropic). 


Volume  fraction,  clftt)  Vokimc  fraction, 


Fig.  3.  Random  SWCNTs  (isotropic).  Fig.  4.  Random  SWCNTs  (isotropic). 
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For  the  nanocomposites  reinforced  by  transversely  isotropic  nanotubes,  the  trends  of  the 
Young’s  modulus  and  Poisson’s  ratio  varying  with  AR  and  c \  (see  Figs.  5-8)  are  similar 
to  those  for  the  nanocomposites  containing  isotropic  nanotubes  (see  Figs.  1  -  4).  An 
exception  is  that  the  Poisson’s  ratio  ( v)  for  isotropic  nanocomposites  reinforced  by 
randomly  oriented  transversely  isotropic  nanotubes  (see  Fig.  8)  decreases  monotonically 
with  the  increase  of  c\,  which  differs  from  what  is  exhibited  by  isotropic  nanocomposites 
reinforced  by  randomly  oriented  isotropic  nanotubes  (see  Fig.  4). 


Fig.  5.  Aligned  SWCNTs  (transversely  isotropic).  Fig.  6.  Aligned  SWCNTs  (transversely  isotropic). 


Fig.  7.  Random  SWCNTs  (transversely  isotropic).  Fig.  8.  Random  SWCNTs  (transversely  isotropic). 


The  stress  transfer  problem  for  traditional  fiber-reinforced  composites  has  been 
extensively  studied.  The  shear-lag  model  originally  proposed  by  Cox  (1952)  provides  a 
good  estimate  of  the  stresses  in  the  fiber  transferred  from  the  matrix  through  the 
interface.  However,  this  analytical  model  cannot  be  directly  applied  to  characterize 
nanotube-reinforced  composites,  since  it  considers  the  load  transfer  across  the  curved 
interface  only  and  regards  the  fiber  ends  as  traction-free.  In  addition,  some  important 
morphological  features  of  nanotube-reinforced  composites  are  not  incorporated  in 
existing  shear-lag  models.  It  is  known  that  the  efficiency  of  the  nanotube  reinforcement 
depends  sensitively  on  the  morphology  (including  diameter,  wall  thickness  and  chirality) 
and  distribution  of  the  nanotubes  (e.g.,  Thostenson  and  Chou,  2003).  This  necessitates  the 
incorporation  of  atomistic  structures  of  nanotubes  in  developing  continuum-based 
analytical  models  of  shear-lag  type  for  nanotube-reinforced  composites.  One  objective  of 
the  current  project  is  to  develop  such  a  shear-lag  model  using  a  representative  volume 
element  (RVE)  of  a  concentric  composite  cylinder  embedded  with  a  capped  carbon 
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nanotube  (see  Fig.  9(b)).  The  molecular  structural  mechanics  of  Li  and  Chou  (2003a)  is 
first  employed  to  predict  the  effective  Young’s  modulus  of  the  capped  nanotube.  The 
capped  nanotube  is  then  replaced  by  an  effective  (solid)  fiber,  whose  length  and  diameter 
are  kept  to  be  the  same  as  those  of  the  nanotube  to  preserve  the  essential  morphological 
features  of  the  nanocomposite  (see  Fig.  9(c)).  The  Young’s  modulus  of  the  effective  fiber 
is  subsequently  determined  under  an  isostrain  condition.  This  is  followed  by  the 
development  of  the  shear-lag  model  using  the  cylindrical  RYE  shown  in  Fig.  9(c). 


Nanotube  (with  radius  a)  Matrix  (with  radius  R) 


(C) 

Fig.  9.  Representative  volume  element. 


The  formulation  of  the  shear-lag  model  is  based  on  the  elasticity  theory  for 
axisymmetric  problems.  The  major  assumptions  are  typical  of  a  shear-lag  analysis.  The 
results  give  (see  Gao  and  Li  (2005)  for  details) 
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By  using  the  shear-lag  formulas  given  in  Eq.  (3),  sample  numerical  results  are  readily 
obtained.  Fig.  10  shows  how  the  normalized  interfacial  shear  stress  and  average  axial 
normal  stress  in  the  nanotube  vary  along  the  nanotube  length  for  three  cases  having 
different  nanotube  aspect  ratios  (ARs).  A  comparison  shows  that  in  all  of  the  three  cases 

(with  different  aspect  ratios)  considered  the  trends  of  both  tJct  and  g[zI  a  varying  with 
z/(2a),  as  predicted  by  the  current  shear-lag  model,  are  the  same  as  those  predicted  by  the 
earlier  computational  model  of  Li  and  Chou  (2003b),  thereby  verifying  the  new  model. 


Em  =  2.41  GPa 
vm  =  0.35 
£>=1000  GPa 
a  =  0.471  urn 
R  =  5a 


Fig.  10.  Interfacial  shear  stress  and  average  axial  normal  stress  in  the  nanotube. 


Fig.  1 1  shows  the  distribution  of  the  average  axial  normal  stress  in  the  nanotube  along 
the  nanotube  length,  while  Fig.  12  illustrates  that  of  the  average  axial  normal  stress  in  the 
matrix.  Clearly,  it  is  observed  from  Fig.  1 1  that  the  maximum  axial  normal  stress  is 
reached  in  the  middle  of  the  nanotube,  whereas  the  minimum  occurs  at  its  two  ends.  Also, 
it  is  seen  that  the  larger  the  aspect  ratio,  the  higher  the  average  axial  normal  stress  at  the 
same  z  for  given  a  and  a.  The  opposite  trends  are  seen  for  the  average  normal  stress  in 
the  matrix,  as  shown  in  Fig.  12.  More  importantly,  it  can  be  noticed  from  comparing  Figs. 
11  and  12  that  when  the  nanotube  aspect  ratio  is  sufficiently  large,  most  of  the  applied 
load  in  the  axial  direction  can  be  taken  up  by  the  nanotube  and  the  surrounding  matrix 
needs  to  support  only  a  small  portion  of  the  axial  loading.  For  example,  when  AR  =  12.8 
the  portion  of  the  applied  load  taken  by  the  nanotube  on  the  cross  section  z  =  0  is  more 
than  three  times  as  large  as  that  by  the  matrix  on  the  same  cross  section,  although  the 
cross-sectional  area  of  the  matrix  is  twenty-four  times  as  big  as  that  of  the  effective  fiber 
on  z  =  0  in  the  current  RVE  with  R  =  5a.  However,  this  is  no  longer  the  case  if  the 
nanotube  aspect  ratio  becomes  small.  For  instance,  when  AR  =  7.8  both  of  the  nanotube 
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and  the  matrix  on  the  cross  section  z  -  0  take  about  the  equal  amount  of  the  applied  load. 
These  reveal  that  the  nanotube  aspect  ratio  is  a  critical  controlling  parameter  for 
nanotube-reinforced  composites,  and  that  for  significant  reinforcements  nanotubes  with 
large  aspect  ratios  should  be  used.  This  observation  conforms  to  what  has  been  found  in 
existing  experimental  and  computational  studies  (e.g.,  Frankland  et  al.,  2003),  thereby 
further  supporting  the  newly  developed  analytical  model. 


Fig.  11.  Average  axial  normal  stress  in  the  nanotube.  Fig.  12.  Average  axial  normal  stress  in  the  matrix. 
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